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APPLICATION  OF  SILICONE  FLUID  RHEOGONIOMETER  MEASUREMENTS 
TO  DESPIN  MOMENT  STUDIES  IN  A  SPINNING  AND  CONING  CYLINDER 


1.  INTRODUCTION 

Liquid  filled  projectile  instad^ility  studies  by  Weber  and  Miller^ 
confirmed  that  high-viscosity  liquid  payloads  can  produce  significant  despin 
moments  that  indicate  the  presence  of  a  destcUailizing  moment  on  the  shell. 
Theoretical  studies^  for  Newtonian  liquids  have  established  that  the  despin 
moment  is  a  function  of  the  liquid  viscosity  and  has  a  maximxim  at  a  low 
Reynolds  number  that  corresponds  to  a  liquid  of  high  viscosity.  The  liquids 
used  by  Weber  cuid  Miller^  in  their  test  fixture  were  Dow  Coming  200  silicone 
fluids  mcuiufactured  by  the  Dow  Chemical  Compeuiy  (Midland,  MI)  .  These  liquids 
are  considered  viscosity  standards  in  many  applications  and  can  be  obtained  in 
specific  viscosities.  They  may  also  be  blended  to  obtain  a  wide  viscosity 
range.  However,  studies  presented  by  other  workers^ have  shown  that  sili¬ 
cone  fluids  exhibit  viscoelastic  behavior.  Weber  and  Miller's  test  results 
for  silicone  fluids^  showed  good  agreement  with  the  theory  developed  for 
Newtonian  fluids  at  a  high  Reynolds  number  (i.e.,  low-viscosity  liquids);  how¬ 
ever,  anomalies  were  seen  at  low  Reynolds  number  for  higher-viscosity  liquids. 

To  determine  the  possible  effect  of  non-Newtonian  liquid  behavior 
of  Dow  Coming  200  silicone  fluids  on  the  despin  moment,  rheology  analysis  was 
performed.  A  Weissenberg  rheogoniometer  was  used  to  measure  the  viscosity  and 
normal  stress  of  the  silicone  fluids  as  a  function  of  shear  rate.  These  data 
were  then  used  to  confute  the  resulting  moments  for  these  fluids  in  a  spinning 
and  coning  cylinder  (simulating  Weber  and  Miller's  experiment^) .  The  results 
indicate  the  presence  cind  influence  of  non-Newtonioui  viscoelastic  effects. 


2  .  WEISSENBERG  RHEOGONIOMETER 

A  Weissenberg  Rheogoniometer,  Model  R18,  was  used  to  measure  the 
rheological  properties  of  the  Dow  Coming  200  silicone  fluids.  Only  a  brief 
description  of  the  Weissenberg  unit  will  be  given  here.  A  more  in-depth 
description  cam  be  found  in  the  instrument's  instructional  manual.^  An 
excellent  review  of  the  Weissenberg  vmit  has  also  been  made  by  Walters^  in 
which  he  discusses  the  instrument  and  its  limitations. 

In  this  study,  the  unit  was  used  in  the  cone  amd  plate  steady 
rotational  configuration  (Figure  1) .  In  this  mode,  stress  measurements  are 
theoretically  possible  for  shear  stress  (<^12)  normal  stress  (N^^)  at 

constamt  shear  rates  (y)  up  to  5, 000/s* However,  in  practice,  the  maiximum 
2shear  rate  at  which  one  cam  obtain  reliadDle  data  is  much  lower  and  dependent 
on  inertia  and  temperature  effects  of  the  liquid  under  test.  Rheology  data 
was  obtained  for  10  silicone  fluids.  The  fluids  were  usually  identified  in 
terms  of  kilo  centistokes  (kcSt)  (e.g.,  600  kcSt  =  600,000  cSt) .  The 
kinematic  viscosities  tested  included  0.1,  1,  4.2,  6,  10,  32,  60,  100,  300, 
amd  600  kcSt  fluids,  encompassing  the  complete  range  of  viscosities  used  in 
the  Weber  amd  Miller  experiments.^ 
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Figure  1 .  Weissenberg  Configuration 


3. 


WEISSENBERG  TEST  RESULTS 


Viscosity  and  normal  stress  test  results  for  the  Dow  Corning  200 
silicone  fluids  are  shown  in  Figures  2  and  3,  respectively.  The  start  of 
shear  thinning  (Figure  2)  is  indicated  at  shear  rates  between  10  and  100/s“^ 
for  higher  viscosity  fluids.  Viscosity  data  for  these  fluids  dropped 
dramatically  (not  shown  in  Figure  2  for  clarity,  but  will  be  shown  later  on 
individual  plots  of  predictive  equations)  at  higher  shear  rates.  This  was 
attributed  to  inertial  or  temperature  effects.  A  normal  stress  was  also 
measured  (Figure  3)  for  fluid  viscosities  as  low  as  4.2  kcSt.  The  normal 
stress  for  higher  viscosity  fluids  also  changed  dramatically  (not  shown  in 
Figure  3  for  same  reason  as  stated  for  Figure  2)  ,  approaching  a  constcint 
stress  at  higher  shear  rates.  Based  on  this  data,  it  was  concluded  that 
non-Newtonian  effects  were  significant  for  higher  viscosity  silicone  fluids, 
which  may  be  responsible  for  the  reported  cuiomalies  in  despin  moments  for 
these  fluids.  To  verify  this,  the  effect  of  viscoelastic  fluid  behavior  on 
despin  moment  was  investigated. 


4  .  VISCOELASTIC  EFFECTS  ON  DESPIN  MOMENT 

Previous  experiments^  and  theoretical  studies®  included  the  effect 
of  viscoelastic  liquid  properties  on  the  despin  moment.  The  theoretical 
studies  made  by  Rosenblat,  Gooding,  and  Engleman®  calculated  the  despin  moment 
for  simple  viscoelastic  models  using  two  approaches.  The  more  rigorous,  but 
more  complex,  technique  involved  calculations  using  a  finite  element  method. 
The  other  technique,  which  was  the  one  used  in  this  study,  involved 
calculations  based  on  assuming  an  infinitely  long  cylinder.  This  method  was 
also  used  by  Herbert^  for  calculating  theoretical  despin  moments  for  Newtonicui 
liquids.  The  calculations  do  not  give  exact  results  for  all  test  conditions. 
However,  good  accuracy  is  obtained  when  applied  to  cylinders  with  small  coning 
to  spin  rate  ratios,  which  was  the  case  in  Weber  and  Miller's  test  fixture.^ 

A  Basic  Language  computer  program  was  written  using  Rosenblat,  Gooding,  and 
Engleman's  derived  equation®  for  the  non-Newtonian  despin  moment.  This 
computer  program  was  then  used  to  determine  the  viscoelastic  effect  of  high- 
viscosity  Dow  Coming  200  silicone  fluids  on  the  despin  moment,  impropriate 
pareuneters  for  simple  viscoelastic  models  used  by  Rosenblat,  Gooding,  and 
Englemam®  were  obtained  from  the  Weissenberg  test  data. 


5.  VISCOELASTIC  MODELS 

Rosenblat,  Gooding,  and  Engleman®  indicated  the  difficulty  in 
determining  the  correct  model  to  use  for  describing  the  non- Newtonian  flow. 
There  is  no  constitutive  relationship  that  can  be  applied  to  all  flows  of  all 
non -Newtonian  liquids.  However,  a  constitutive  equation  describing  specific 
flows  may  be  applicadDle  for  a  wide  range  of  fluids.  Rosenblat,  Gooding,  and 
Engleman®  also  indicated  that  ”it  is  not  known  whether  shear  thinning,  noirmal 
stress  differences,  elongational  viscosity,  stress  relaxation,  or  some 
combination  of  all  or  some  of  these”  is  significant  in  any  real-flow  problem. 
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Figure  2.  Dow  Coming  200  Viscosity  Results 
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LOG  1/SECONDS 


Figure  3.  Dow  Corning  200  Normal  Stress  Results 


11 


Therefore,  Rosenblat,  Gooding,  and  Engleman®  restricted  themselves  to  testing 
and  comparing  simple  constitutive  relationships  to  experimental  data. 

Rosenblat,  Gooding,  and  Engleman^  selected  two  viscoelastic  models 
to  test;  a  "second  order  fluid”  and  a  "differential  model.”  In  a  second  order 
fluid,  the  following  relationships  are  valid: 


n(Y)  =  no 

Ni(y)  =  2X»ioY^ 


i1q  is  the  zero  shear  viscosity,  cuid  1  is  a  relaxation  time  constant 
for  the  fluid.  As  seen,  shear  thinning  is  not  predicted,  and  normal  stress  is 
a  quadratic  fxinction  of  shear  rate.  Weissenberg  data  for  higher-viscosity, 
silicone  fluids  at  high-shear  rates  did  indicate  shear  thinning;  therefore, 
this  model  does  not  appear  to  be  applicable. 

For  the  differential  model,  a  ”3  constcuit  modified"  version  of  the 
"Oldroyd  8  constant  model was  used.  The  3  constant  modified  model  exhibits 
both  shear  thinning  auid  normal  stress  behavior.  The  3  constant  modified  model 
includes  a  zero- shear  viscosity  tiq,  a  relaxation- time  constant  1,  and  a 

retardation- time  constant  where  e  defines  the  hi  ratio. 

X 


The  pertinent  equations  are: 


(l-€) 


A  Basic  Language  computer  program  was  written  (Appendix  A)  to 
calculate  the  fluid  relaxation  time  for  four  of  the  Dow  Coming  200  silicone 
fluids.  The  four  silicone  fluids  corresponded  to  fluids  used  extensively  by 
Weber  and  Miller^  and  included  600,  300,  100,  cuid  10  kcSt  fluids.  The 
relcucation-time  constant  for  each  fluid  was  determined  using  the  Weissenberg- 
zero  shear  viscosity  cuid  normal  stress  versus  shear  rate  data.  The  zero- shear 
viscosity  is  calculated  directly  from  the  Weissenberg  measurements.  Normal 
stress  data  at  the  lower  shear  rates  were  substituted  into  the  normal  stress 
equation,  which  was  then  solved  for  X  at  assumed  values  of  €  =  0,  .2,  and  .4. 
An  average  value  for  X  was  then  obtained  from  these  results.  An  exact 
determination  of  e  for  silicone  fluids  could  not  be  made  because  of 
limitations  in  Weissenberg  measurements  at  the  high-shear  rates.  Therefore, 
assumed  values  of  €  were  used  to  determine  fluid  relaixation  times.  Table  1 
shows  the  results  of  these  calculations. 
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Table  1.  Calculated  Silicone  Fluid  Relaxation  Times 


Viscosity  Relaxation  Time 

kcSt  6=0  €=.2  €=.4 


600 

.0144 

.018 

.0244 

300 

.0099 

.0124 

.0166 

100 

.003 

.0038 

.0051 

10 

.00138 

.00173 

.00233 

Predictions  of  viscosity  and  normal  stress  obtained  for  the  3  con¬ 
stant  modified  model  as  a  function  of  shear  rate  are  shown  in  Figures  4-7  for 
e  =  0.  (NOTE:  The  viscosity  prediction  for  the  other  values  of  €  indicates  a 
constant  viscosity  of  lower  magnitude  after  a  short- shear,  thinning  interval. 
The  predicted  values  of  normal  stress  also  become  constant  after  this 
interval.)  All  test  data  points  are  included  in  the  plots;  a  good  data  fit  is 
indicated.  These  fluid  constants  were  then  used  to  calculate  the  despin 
moment  using  the  equation  developed  by  Rosenblat,  Gooding,  and  Engleman®  for 
the  differential  model. 


6.  DESPIN  ^MENT  CALCULATIONS 

The  despin  moment  equation  is  shown  below: 

M  =  ■  ■— — -  =  REAL  I  1 


where 


M  =  normalized  despin  moment 
3=  actual  despin  moment 
A  =  coning  rate  to  spin  rate  ratio 
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Figure  4 .  10  kcSt  Rheology  Predictions 
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100  KCS  SILICONE  FLUID 
PREDICTED  VISCOSITY  VERSUS  SHEAR  RATE 
4  -TEST  DATA  POINTS 
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Figure  5.  100  kcSt  Rheology  Predictions 
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300  KCS  SILICONE  FLUID 
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Figure  6.  300  kcSt  Rheology  Predictions 
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♦  -TEST  DATA  POINTS 


-10  1  2  3  4 


LOG  1/SECOHOS 


600  KCS  SILICONE  FLUID 
PREDICTED  NORMAL  STRESS  VERSUS  SHEAR  RATE 
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Figure  7 .  600  kcSt  Rheology  Predictions 
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B  =s  cylinder  aspect  ratio  (4.5  for  the  reported  tests) 


C  =  sine  of  the  coning  angle 

The  term  in  brackets  is  a  solution  of  the  differential  equation  for  the  fluid 
flow  field,  which  is  a  function  of  the  fluid  model  specified  by  the  variable 
S.  The  varicdDle  S  for  the  various  models  is  shown  below: 

S=Re  --  Newtonicui  model 

S=Re/ (1+iDeRe)  --  second  order  fluid  model 
S= (1-iDeRe) Re/ (1-  iDeRe)  --  differential  model 


Where  Re  is  a  Reynolds  nxjmber  curid  De  is  a  Deborcih  niomber,  both  dimensionless 
parameters  defined  by  Rosenblat,  Gooding,  and  Engleman.^  Reynolds  number  is 
proportional  to  the  spin  rate  of  the  cylinder  and  inversely  proportional  to 
the  fluid  zero-shear  viscosity.  The  Deborah  number  is  proportional  to  the 
product  of  the  fluid  relaxation  time  and  zero- shear  viscosity. 

A  Basic  Language  computer  program  (Appendix  B)  was  written  using 
the  despin  moment  equation.  The  software  was  written  in  a  user-friendly 
format  that  required  information  be  entered  as  input  data  on  request  by  the 
program.  The  first  input  request  is  the  type  of  fluid  model  (i.e.,  Newtonian, 
second  order,  or  differential) .  The  fluid  name  and  Deborah  number  are  then 
requested  for  viscoelastic  fluids  in  the  next  two  user  inputs.  Table  2  shows 
the  calculated  Deboradx  ntimbers  for  the  various  silicone  fluids. 

Tahle  2.  Calculated  Silicone  Fluid  Deborah  Number 


Viscosity 

Deboradi  Nxjmber 

kcSt 

6  =  0 

6  =  .2 

6  =  .4 

600 

2.96 

3.70 

5.01 

300 

1.01 

1.26 

1.69 

100 

0.105 

0.133 

0.179 

10 

0.0048 

0.0061 

0.008 

If  a  differential  fluid  model  is  being  used,  the  program  then  asks 
for  the  €  value  for  the  model  equation.  At  this  point,  all  required  fluid 
specifics  have  been  entered  into  the  program. 
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The  next  information  requested  is  the  type  of  despin  moment  to  be 
calculated  (i.e.,  linear  or  generalized).  This  requires  simple  multiplication 
by  an  appropriate  factor.  For  a  linearized  moment,  the  factor  is  equal  to 
four  times  the  coning  rate  to  the  spin  rate  ratio.  For  the  generalized 
moment,  the  factor  is  equal  to  four  times  the  square  of  the  cosine  of  the 
coning  angle.  The  Reynolds  nxamber  range  and  increments  for  making  the 
calculations  are  then  requested.  These  requests  were  included  in  the  program 
for  two  reasons.  First,  the  i^ple  lie  ccxnputer  used  in  this  study  has 
insufficient  memory  (64K)  to  perform  calculations  over  the  full  rcuige  of 
Reynolds  numbers  used  by  Rosenblat,  Gooding,  and  Englemaui.®  Second,  and  more 
significant,  the  calculated  values  of  despin  moment  would  cover  only  the  range 
of  Reynolds  ntimber  used  in  the  test  program.  The  test  data  determined  by 
Weber  atnd  Miller^  was  obtained  over  a  limited  ramge  of  Reynolds  number  because 
of  limitations  in  the  spin  rate  of  their  test  fixture.  The  remaining  program 
requests  information  pertaining  to  obtaining  copies  of  the  calculated  despin 
moments . 


7.  DISCUSSION  OF  DESPIN  MOMENT  RESULTS 

Figure  8  shows  the  test  results  obtained  by  Weber  and  Miller^  for 
various  silicone  fluids.  Theoretical  curves  obtained  by  Herbert^  and 
Rosenblat,  Gooding,  and  Engleman®  for  the  despin  moments  are  include  in  Figure 

8.  The  data  points  that  deviate  from  the  theoretical  curve  at  the  low 
Reynolds  numbers  correspond  to  silicone  viscosities  of  600,  300,  and  100  kcSt, 
from  left  to  right.  If  curves  are  drawn  through  the  data  points  for  each 
viscosity,  a  family  of  curves  can  be  established  for  these  three  silicone 
viscosities  at  the  low  Reynolds  numbers.  The  10  kcSt  data  points,  as  well  as 
other  lower  silicone  viscosity  fluids  tested,  appear  to  be  close  to  the 
theoretical  prediction  of  despin  moment. 

Figures  9,  10,  and  11  show  the  predicted  despin  moments  using  the 
differential  model  parameters  shown  in  Tables  1  and  2 .  The  three  curves 
correspond  to  e  equal  to  0,  0.2,  and  0.4.  Included  in  each  Figure  is  the 
despin  moment  response  for  a  Newtonian  fluid.  (NOTE:  The  response  for  the 
Newtonian  fluid  indicates  a  maximum  at  approximately  a  Reynolds  nxjmber  of  15 
con^ared  to  a  Reynolds  number  of  approximately  40  in  Figure  8  for  the 
theoretical  curves  by  Herbert^  and  Rosenblat,  Gooding,  cind  Engleman,^  This 
difference  is  due  to  the  infinite  cylinder  approximation  used  in  this 
auialysis.  However,  the  general  shapes  of  the  curves  are  the  same). 

For  e  =  0  (Figure  9) ,  the  despin  moment  response  for  the  600  kcSt 
fluid  rises  sharply  near  a  Reynolds  niimber  of  unity  and  returns  to  an 
approximate  value  of  0.2  at  Reynolds  number  3  (this  rauige  of  1-3  corresponds 
to  the  actual  Reynolds  number  range  used  in  the  study) .  The  300  kcSt  fluid 
has  two  peaks  in  the  2-8  Reynolds  numbers  ramge;  whereas,  the  100  kcSt  fluid 
has  one  peak  in  the  5-20  Reynolds  numbers  ramge.  In  all  three  cases,  the 
peaks  extend  well  beyond  the  value  of  despin  moment  obtained  for  a  Newtoniam 
fluid.  The  10  kcSt  fluid  results  appear  somewhat  higher  tham  the  Newtonian 
prediction.  However,  it  is  decreasing  as  e3q)ected  in  this  Reynolds  nximber 
ramge  (50-100) .  These  results  agree  with  the  trend  in  the  experimental  data, 
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which  show  significantly  larger  despin  moments  for  the  higher  viscosity 
fluids;  whereas,  Newtonian  theory  predicts  much  lower  despin  moments. 


For  the  other  two  values  of  €,  €  =  0.2  (Figure  10)  and  €  =  0.4 
(Figure  11) ,  the  peaks  mentioned  in  the  above  paragraph  are  significantly 
reduced.  The  resultant  curve  for  the  100  kcSt  fluid  qualitatively  approaches 
that  of  the  test  data,  whereas,  the  curves  for  the  300  and  600  kcSt  exhibit 
lower  despin  moments  than  the  experimental  test  data.  However,  in  all  cases, 
the  differential  model  predictions  are  higher  than  the  despin  moments  for  a 
Newtonian  fluid. 

It  appears  that,  based  on  the  above  cuialysis,  a  differential  model 
can  be  used  to  predict  the  qualitative  differences  in  despin  moment  for  the 
high- viscosity,  Dow  Coming  200  silicone  fluids  at  a  low  Reynolds  number. 
However,  precise  quantitative  agreement  with  test  data  is  not  indicated. 

Much  better  agreement  with  test  data  can  be  found  if  one 
substitutes  a  Deborah  number  ec[ual  to  0.3  times  the  calculated  values  obtained 
for  €  s  0 .  Figure  12  shows  the  predicted  despin  moments  for  this  condition. 
Comparison  of  the  despin  moments  for  the  300  and  600  kcSt  fluids  with  Weber 
and  Miller's  test  data^  (Figure  8)  shows  a  sharp  increase  with  increasing  the 
Reynolds  number.  The  100  kcSt  fluid  shows  a  rise  and  peak  at  a  Reynolds 
number  of  edDout  20.  The  10  kcSt  shows  a  decreasing  trend  that  is  very  close 
to  the  Newtonicui  prediction.  All  these  fluids  behave  similarly  to  the  actual 
test  results. 


8 .  CONCLUSIONS 

A  Weissenberg  rheogoniometer  was  used  to  measure  the  rheological 
properties  for  a  series  of  Dow  Coming  200  silicone  fluids.  The  following 
information  was  foimd  from  these  measurements: 

•  Non-Newtonicui  behavior  was  indicated  for  the  higher  viscosity 
silicone  fluids  in  the  form  of  shear  thinning  and  normal  stress  at  high  shear 
rates . 


•  "Differential  model”  viscoelastic  parameters  could  be 
determined  from  rheogoniometer  measurements.  A  "second  order  fluid  model”  was 
found  to  be  inappropriate  because  fluid  shear -thinning  effects  are  not 
predicted  by  this  model . 

•  The  "differential  model"  qualitatively  predicted  the  despin 
moment  behavior  of  the  silicone  fluids.  However,  quantitative  agreement  was 
not  achieved. 

•  Much  better  agreement  with  test  data  was  fotind  for  all  fluids 
if  the  "differential  model”  viscoelastic  parameter  for  Deboredi  number  was  made 
equal  to  0.3  times  the  calculated  value  obtained  for  £  :=  0 . 
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Figure  8.  Weber  and  Miller  Despin  Moment  Results 
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Figure  9.  Predicted  Despin  Moment  {e  =  0) 
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Figure  10.  Predicted  Despin  Moment  (6  =  .2) 
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Figure  11.  Predicted  Oespin  Moment  (e  s  .4) 
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APPENDIX  A 


FLUID  RELAXATION  TIME  -  BASIC  LANGUAGE  PROGRAM 


0 

TEXT  :  HOME 

2 

ONERR  GOTO  77 

5 

INPUT  “VISCOSITY" ;N0 

7 

INPUT  “N1,R1  "jNl,Rl 

3 

INPUT  ''N2,R2  ";N2,R2 

10 

INPUT  "EhITER  E  VALUE  ";E 

1  1 

Z  =  2  -  <2  *  E) 

12 

INPUT  "R,S,T  "  ;R,S,T 

15 

FOR  K  =  R  TO  8  STEP  T 

16 

A1  =  1 :A2  =  1 

18 

B1  =  -  <  <Z  *  NO)  /  <N1  * 

K)  ) 

19 

Cl  =  1  /  <K  »  R1  *  Rl) 

20 

Y1  =  -  B1  +  <<B1  *  Bl)  - 

A1  *  CD)  *  .5 

<4  * 

30 

Y1  =  Y1  /  <2  *  Al) 

40 

XI  =  -  Bl  -  (<B1  *  Bl)  - 

Al  »  Cl ) )  *  .5 

<4  * 

50 

XI  =  XI  /  (Al  *  2) 

52 

B2  =  -  < (Z  *  NO)  /  <N2  » 

K)  ) 

53 

C2  =  1  /  <K  *  R2  »  R2) 

60 

Y2  =  -  B2  +  <<B2  *  B2)  - 

A2  *  C2))  *  .5 

(  4  * 

64 

Y2  =  Y2  /  (2  »  A2) 

66 

X2  =  -  B2  -  <(B2  »  B2)  - 

A2  *  C2))  *  .5 

<4  » 

68 

X2  =  X2  /  <A2  *  2) 

69 

PRINT  K 

70 

PRINT  Y1 ,X1 :  PRINT  Y2,X2 

72 

PRINT 

73 

INPi.iC  A* 

75 

NEXT  K 

77 

INPUT  "HIT  KE'r  FOR  NEXT  K 
GE  " ;A« 

RAN 

30 

GOTO  12 

29 
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APPENDIX  B 


DESPIN  MOMENT  -  BASIC  LANGOAGE  PROGRAM 


1  V  =  0  ;  K  =  0  ;  PU)  =  .5 
5  DIM  T1 < 100) 

10  TEXT  :  HOME  ;  (vTAB  <5);  PRINT 
"SELECT  FLUID  MODEL  (1,2  OR 
3)  " 

1 2  MTAB  (10) 

13  PRINT  TAB(  6>"1.  NEWTONIAN  F 

LUID":  PRINT 

14  PRINT  TAB(  6) "2.  SECOND  ORDE 

R  FLUID":  PRINT 

15  PRINT  TAB(  6) "3.  DIFFERENTIA 

L  MODEL  FLUID":  PRINT  :  PRINT 
;  PRINT 

16  PRINT  :  INPUT  AZ« 

17  IF  A2^  =  "1"  THEN  DE%  =  "0":E 

$  =  "0":A2^  =  "NEWTONIAN  FLU 
ID":DE  =  UAL  <DE^):E  =  UAL 
(E*):NF$  =  "NEWTONIAN":  GOTO 
50 

18  IF  A2^  =  "2"  THEN  £%  =  "0":PW 

=  -  PW:E  =  UAL  (E$);  GOTO 

30 

19  IF  A2S  =  "3"  THEN  GOTO  30 

20  GOTO  1 6 

30  PRINT  :  INPUT  "ENTER  FLUID  NA 

ME  " ;NFS 

31  PRINT  :  INPUT  "ENTER  DEBORAH 

NUMBER  ";DE*;DE  =  UAL  (DE*) 

32  IF  AZ^  =  *2"  THEN  AZ*  =  “SECO 

ND  ORDER  FLUID";  GOTO  50 
34  PRINT  ;  INPUT  "ENTER  E  UALUE 
" ;E^:E  =  UAL  <E«) 

36  AZ^  =  "DIFFERENTIAL  MODEL  FLUI 
D" 

50  PRINT  ;  INPUT  "LINEAR  OR  GENE 
RAL  DESPIN  MOMENT  <L  OR  G>  ? 

"  ;AX* 

52  IF  AXS  =  "L"  THEN  AY$  =  "LINE 
ARIZED";  GOTO  60 

54  IF  AX»  =  "G"  I  HEN  AYS  =  "GENE 
RAL I  ZED";  GOTO  60 
56  GOTO  50 

60  PRINT  ;  INPUT  "ENTER  LOWER  RE 
YNOLDS  NUMBER  " ; SLS ; SL  *  UAL 
(SLS) 

65  PRINT  ;  INPUT  "ENTER  UPPER  RE 
YNOLDS  NUMBER  " ; SUS ; SU  =  UAL 
(SUS) 

70  PRINT  ;  INPUT  "ENTER  REYNOLDS 
NUMBER  INCREMENT  ";SIS;SI  = 
UAL  (SIS) 
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72  PRINT  !  INPUT  "DO  YOU  WANT  A 

HARDCOPY  OF  THE  DATA?  CY  OR 
N)  "  iiA$ 

73  IF  A«  =  "Y“  THEN  PRINT  CHR« 

<4) ; "PR#1 " :  GOTO  75 

74  IF  A$  <  >  "N“  THEN  GOTO  72 

75  HOME 

SO  C$  =  AZt;  60SUB  600 
81  IF  DE  <  >0  THEN  =  NF^ :  GOSUB 

600 

92  IF  DE  <  >0  AND  LEN  (A2#)  = 

24  THEN  PRINT  TABC  9)"DEB0 
rah  #=  “DE"  E=  "E 
83  IF  DE  <  >0  AND  LEN  (AZ^)  = 

IS  THEN  PRINT  TAB(  14) "DEB 
ORAH  #=  "DE 

85  C*  *  AY^:  GOSUB  600 

87  IF  AX^  =  "G"  THEN  PRINT  TAB< 

12) "DESPIN  MOMENT  FOR"!  PRINT 
TAB(  12) "CONING  ANGLE=" :  GOTO 
90 

89  PRINT  TAB<  12) "DESPIN  MOMENT 

FOR  T=" 

90  PRINT  TAB(  0 J " - 


91  IF  AXS  =  "G"  THEN  GOSUB  500 

92  IF  AX«  =  "G"  THEN  GOTO  96 

93  PRINT  "RE  ;  TAB<  8) "0.2";  TAB< 

16) "0.1 5";  TAB<  24) "0.1";  TAB< 
32) "0.05" 

95  PRINT  " - ";  TAB<  8)" - ";  TAB< 

16)" - ";  TAB<  24)" - ";  TAB< 

32) " - " 

96  PRINT 

99  WP  *  0:RN  =  0 : IN  =  0:RD  =  0:ID 

*  0 

100  FOR  S  =  SL  TO  SU  STEP  SI 

105  *  V  +  1 

120  P  a  LOG  <<2  »  K)  +  4)  /  LOG 
<10):R  =  (<2  #  K)  +  1)  »  LOG 
<2)  /  LOG  <10) 

1 22  X  a  < < DE  »  DE  *  S  *  S)  +  1)  * 

PW 

123  X  a  s  #  X 

125  Ya  <<DE*DE*S«S*E*E 
)  +  1)  ‘  .5 

127  Q  a  K  «  LOG  (X  /  Y)  /  LOG  < 

10) 

130  IF  K  a  0  THEN  0  =  0:U)W  =  0:  GOTO 
170 

140  0  a  K  *  <  -  1.57-  ATN  <  DE  * 

S)  +  ATN  (E  #  DE  *  S)) 
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150  WU)  =  WP  +  (  LOG  <K)  /  LOG  <1 
0>  ) 

170  ZZ  =  UIW  +  <  LOG  <K  +  1>  /  LOG 

<  10>  ) 

175  NN  =  Q  -  <P  +  R  +  UIW  +  ZZ>:DD 

=  Q  -  (R  +  WW  +  ZZ) 

176  N  =  10  *  NN;D  =  10  *  DD 

185  RR  =  N  »  COS  (0):II  =  N  *  SIN 

<0) 

187  RS  =  D  *  COS  (0);IS  =  D  »  SIN 

<0) 

210  RN  =  RN  +  RRrIN  =  IN  +  II 

220  RD  =  RD  +  RS;ID  =  ID  +  IS 

221  IF  RN  =  R1  AND  RD  =  R2  THEN 
K  =  K  -  1 :  GOTO  240 

222  IF  K  <  4  THEN  GOTO  224 

223  IF  RN  =  R1  AND  RD  =  R2  THEN 
K  =  K  -  1 :  GOTO  240 

224  R1  =  RN:R2  =  RDtll  =  IN:I2  = 

ID 

225  K  =  K  +  1 
227  WP  =  WW 
230  GOTO  1 20 

240  SN  =  ((RN  »  RN)  +  (IN  *  IN))  * 

.5 

245  IF  RN  <  0  AND  IN  <  0  THEN  00 
=  -3.14+  ATN  ( IN  /  RN) : 

GOTO  260 

247  IF  RN  <  0  AND  IN  >  0  THEN  RN 

=  -  RN:00  =  3.14  -  ATN  (I 

N  /  RN)!  GOTO  260 
250  00  =  ATN  (IN  /  RN) 

260  SD  =  ((RD  *  RD)  +  (ID  *  ID))  * 
.5 

265  IF  RD  <  0  AND  ID  <  0  THEN  OD 
=  -3.14+  ATN  (ID  /  RD) : 

GOTO  280 

267  IF  RD  <  0  AND  ID  >  0  THEN  RD 

=  -  RDrOD  =  3.14  -  ATN  (I 

D  /  RD):  GOTO  280 


270 

OD  =  ATN 

(ID  /  RD) 

280 

ST  =  SN  / 

SD:OT  =  (00  -  OD) 

290 

IM  =  ST  » 

SIN  (OT) 

300 

RE  =  1  » 

IM 

310 

REM  (- 

J*J=1 ) 

315 

IF  AXS  = 

"G"  THEN 

GOSUB  550 

317 

IF  AX*  = 

"G"  THEN 

GOTO  325 
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320 


PRINT  S?  TAB(  8)  INT  <8000  * 

RE)  /  10000;  TAB<  16)  INT  <6 
000  *  RE)  /  10000;  TAB<  24)  INT 
<4000  *  RE)  /  10000;  TAB<  32 
)  INT  <2000  *  RE)  /  10000 
325  T1<M)  =  RE 

330  WP  =  0;RN  =  0 ; IN  =  0 ; RD  =  0;I 
D  =  0 
335  K  =  0 
340  NEXT  S 

345  IF  A^  =  "Y"  THEN  PRINT  CHR^ 
<4) ; "RRMO" 

350  INPUT  "DO  YOU  WANT  TO  SAOE  D 
ATA  TO  DISK  <Y  OR  N)?  “ ;A^ 

355  IF  A#  =  "N"  THEN  END 
360  IF  A»  <  >  "Y"  THEN  GOTO  35 

0 

362  PRINT  ;  INPUT  "ENTER  DRWE  N 
UMBER  FOR  DATA  DISK";DN«;DN  = 

‘v)AL  <DN«) 

363  IF  DN  =  2  OR  DN  =  1  THEN  GOTO 
365 

364  GOTO  362 

365  IF  DE  =  0  AND  E  =  0  THEN  = 

NFS  +  +  SLS  +  +  BUS  + 

+  SIS:  GOTO  380 

367  IF  MIDS  <AZS,1,1)  =  "S"  THEN 
FS  =  NFS  +  +  DES  +  + 

SLS  +  ^  SUS  +  +  SIS; 

E  =  -  1 ;  GOTO  380 

370  FS  =  NFS  +  +  dES  +  + 

ES  +  +  SLS  +  +  SUS  + 

+  SIS 

380  POKE  43624, DN 
385  DS  =  CHRS  <4) 

390  PRIN)  i>S;"OPEN  "FS 

400  PRINT  DS; "WRITE  "FS 

405  PRINT  NFS 

410  PRINT  U 

412  PRINT  DE 

414  PRINT  E 

416  PRINT  SL 

420  PRINT  SU 

430  PRINT  SI 

440  FOR  K  a  1  TO  U 

450  PRINT  T1<K) 

470  NEXT  K 

480  PRINT  DS;" CLOSE  "FS 
485  POKE  43624,1 
490  END 
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500  PRINT  "RE  ♦»"  ;  TAB<  8)  "20";  TAB< 
16) “15";  TAB(  24) "10" 

510  PRINT  “ - ";  TAB<  8>" - ";  TAB( 

16)" - ";  TAB<  24)“ - " 

520  RETURN 

550  PRINT  S;  TAB<  8)  INT  (35300  * 

RE)  /  10000;  TAB<  16)  INT  (3 
7300  *  RE)  /  10000;  TAB(  24) 

INT  (38800  «  RE)  /  10000 
590  RETURN 
600  XX  =  LEN  (CS) 

610  YY  =  (40  -  XX)  /  2;YY  =  INT 
(YY) 

620  PRINT  TAB(  YY)C$ 

630  RETURN 
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